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investigated.  If  the  current  force  overcomes  the  sea  bottom  resistance,  the  hose 
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evaluated  in  terms  of  lateral  deflections,  hose  tensions,  and  anchor  loads.  The 
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of  the  simulation  model  are  briefly  summarized.  A  parametric  analysis  is  conducted  to 
study  the  influence  on  the  hose  response  of  the  physical  parameters  considered  in  the 
simulation  model.  The  results  indicate  that,  for  a  practical  hoseline,  the  most  critical 
parameters  are:  the  segment  length-to-span  ratio,  the  axial  rigidity  of  the  hose,  the 
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INTRODUCTION 


Objective 

The  objective  of  this  study  was  to  determine  analytically  the  lateral 
stability  of  the  Navy's  flexible  submarine  hoselines  in  a  slowly  varying 
current.  This  goal  was  achieved  through  the  development  of  three  easy 
to  use  design  tools:  a  numerical  simulation  model,  design  charts,  and  a 
parametric  model. 

Scope 

This  effort  analyzed  the  lateral  stability  of  a  flexible  submarine 
hoseline  in  varying  currents.  The  analysis  was  limited  to  a  long  hose 
of  little  flexural  and  torsional  stiffness.  Further,  only  the  behavior 
of  a  taut  hose  segment  was  studied.  The  response  of  a  slack  hose  segment, 
which  has  zero  axial  tension,  was  not  included.  Theories  for  simulating 
the  behaviors  of  a  submarine  hosel.ine  in  currents  were  derived.  Design 
tools  were  developed  to  provide  the  required  stability  for  engineering 
applications. 

This  report  presents  the  development  of  the  hoseline  simulation 
model.  Pertinent  theories  used  in  the  derivation  of  the  model  are  briefly 
described.  The  numerical  solution  techniques  incorporated  in  the  computer 
code  are  included  in  Appendix  A.  The  procedures  and  major  findings  obtained 
from  the  parametric  analysis  are  also  discussed.  Results  of  the  analysis 
are  presented  in  graphs  to  show  the  influence  of  each  governing  parameter. 
Design  charts  and  parametric  models  for  typical  Navy  hoses  are  presented 
followed  by  a  guideline  for  using  these  results  for  design. 


Background 


Large  volumes  of  fuels  are  consumed  daily  by  amphibious  forces  when 
engaged  on  a  battle  field.  Liquid  cargo  of  the  required  volume  can  be 
transferred  only  by  tankers,  which  may  be  moored  several  miles  from  shore. 
The  cargo  is  discharged  through  conduits  to  the  beachhead.  Time  is  a 
critical  factor;  it  dominates  system  design  in  addition  to  functional 
reliability.  Thus,  the  most  effective  system  is  the  one  that  is  trans¬ 
ported  easily  and  rapidly,  requires  the  least  resources  for  installation 
and  operation,  and  is  largely  universal.  The  Naval  Civil  Engineering 
Laboratory  (NCEL)  is  currently  developing  a  flexible  high-strength  con¬ 
duit  to  maximize  the  effectiveness  of  the  Navy's  liquid  cargo  handling 
system.  The  conduit,  which  is  highly  collapsible  and  allows  a  very 
small  bending  radius,  could  be  compactly  spooled  on  hose  reels  for  easy 
storage  and  transportation. 

Among  the  engineering  considerations,  the  most  important  factor  in 
a  submarine  flowline  design  is  to  ensure  lateral  stability  against  envi¬ 
ronmental  loads.  The  most  frequently  used  stabilization  techniques  in 
industrial  practice  are  weight  coating,  trenching,  and  anchoring.  Heavy 
coating  negates  the  collapsibility  of  the  conduit,  and  commercial  pipe¬ 
line  trenching  operations  require  heavy  equipment  and  intensive  labor  at 
the  site.  These  features  are  obviously  contrary  to  the  primary  objective 
of  expediting  the  military  contingency  operation.  Therefore, 
securing  the  flowline  with  mechanical  anchors  seems  to  be  the  most 
practical  approach  for  stabilization.  There  have  been  a  number  of 
extensive  investigations  on  the  stability  of  semirigid  submarine 
pipeline  systems  (Ref  1,  2,  and  3).  However,  highly  flexible  hoselines 
have  received  little  attention.  A  direct  numerical  integration 
technique  has  been  employed  successfully  to  describe  the  catenary  shape 
of  a  long  cable  subject  to  steady  ocean  current  (Ref  4,  5,  and  6).  This 
effort  was  to  investigate  the  feasibility  of  using  the  same  technique  to 
determine  the  behavior  of  a  cable-like  hoseline  system. 
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SIMULATION  MODEL 


Problem  Definition 

A  definition  sketch  of  a  multisegment  hoseline  on  the  ocean  bottom 
is  given  in  Figure  1.  The  hoseline  may  be  restrained  at  arbitrary  loca¬ 
tions  to  the  bottom  with  mechanical  anchors.  A  pertinent  free-body  dia¬ 
gram  is  given  in  Figure  2.  A  right-handed  coordinate  system  X^,  X2,  X^ 
with  unit  orthogonal  base  vectors  e^,  e^  embedded  in  the  bottom 
with  e^  directed  away  from  the  ocean  bottom.  The  vector  e^  is  directed 
outward  perpendicular  to  the  shore.  The  bottom  is  gently  sloping  in  the 
e^  direction  at  an  angle  </>  with  respect  to  the  ocean  surface.  The  gravity 
vector  k  makes  an  angle  with  e^. 

It  is  assumed  that  the  velocity  vector  on  the  bottom  is  oriented  in 
the  direction  of  There  is  no  e^  component  of  the  velocity  vector, 

and  its  magnitude  may  vary  in  a  piecewise  linear  fashion  with  X^.  The 
velocity  is  assumed  to  be  slowly  varying  (i.e.,  at  any  instant  a  steady 
current  will  be  assumed).  No  inertial  forces  will  be  considered.  The 
velocity  induces  drag  and  lift  forces  on  the  hoseline.  Morison-type 
drag  forces  are  assumed  and  the  independence  principle  is  invoked  to 
separately  specify  drag  coefficients  in  the  normal  and  tangential  direc¬ 
tions  by  multiplying  the  squares  of  the  normal  and  tangent  components  of 
velocity,  respectively. 

The  hoseline  may  comprise  several  segments  joined  end-to-end  and 
forming  a  curved  line  on  the  bottom,  with  each  segment  having  different 
cross-sectional  properties.  There  are  anchor  points  at  the  junctions  of 
the  segments.  It  is  assumed  that  each  anchor  behaves  as  a  linearly  elastic 
spring  acting  in  the  direction  of  the  force  necessary  to  equalize  the 
end-point  tensions  on  the  incident  segments. 

It  is  assumed  that  each  segment  slides  on  the  bottom  without  roll¬ 
ing  or  lifting  off  the  bottom.  An  analysis  including  the  couj.  Ling  between 
axial  tension  and  torsional  responses  was  beyond  the  scope  of  this  work. 

In  fact,  the  torsional  stress,  which  is  much  smaller  than  that  of  axial 
tension,  will  not  significantly  affect  the  catenary  shape  of  the  hoseline 
(Ref  6).  Therefore,  the  current -induced  loads  on  the  hose  are  insensitive 
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to  the  degree  of  twisting,  as  long  as  the  hose  cross  section  remains, 
circular  without  serious  buckling.  Consequently,  the  axial  and  the  tor¬ 
sional  responses  of  the  hoseline  are  assumed  separable  in  this  develop¬ 
ment.  There  is  a  lift  force  from  the  bottom  current  which  tends  to  lift 
the  segment,  but  if  the  segment  were  to  separate  from  the  bottom,  the 
lift  coefficient  would  decrease  and  the  segment  would  be  forced  back 
into  contact  with  the  bottom,  as  discussed  in  Reference  1. 

The  frictional  force  resisting  sliding  of  each  segment  is  assumed 
to  be  of  the  Coulomb  type  and  to  be  separable  into  components  tangent 
and  normal  to  the  segment,  with  each  component  having  a  different 
coefficient  of  friction.  The  friction  resistance  is  directly  pro¬ 
portional  to  the  magnitude  of  the  net  force  normal  to  the  seafloor  by 
the  coefficient  of  static  friction  if  the  segment  is  immobile,  and  by 
the  coefficient  of  kinetic  friction  as  soon  as  the  segment  moves. 

Typical  value  of  the  kinetic  coefficient  is  approximately  25  percent 
smaller  than  that  of  the  static  coefficient. 

Scenario  of  Motion 

With  the  hoseline  resting  on  the  bottom,  equilibrium  can  be 
established  under  zero  velocity  conditions.  Then,  as  the  velocity 
increases,  a  scenario  of  motion  can  be  constructed  as  follows:  (1)  as 
the  velocity  increases,  the  drag  and  lift  forces  on  the  segments  in¬ 
crease,  (2)  the  magnitude  of  the  reactive  force  from  the  bottom  de¬ 
creases  and,  thus,  the  static  Coulomb  friction  decreases,  (3)  the 
friction  force  equilibrating  the  fluid  drag  forces  and  tensions  in  the 
curved  segments  increases,  (4)  if  the  friction  forces  in  a  particular 
segment  exceed  the  holding  capacity  over  a  significant  fraction  of  the 
segment  span,  the  segment  will  move  to  reestablish  an  equilibrium  posi¬ 
tion  under  those  velocity  conditions,  (5)  as  the  segment  moves,  the 
frictional  force  decreases  in  magnitude  to  that  predicted  using  dynamic 
friction  coefficients  and  has  friction  opposing  the  motion  (i.e,  in  the 
same  direction  as  had  the  static  frictional  resistance),  (6)  new  loca¬ 
tions,  orientations,  and  tension  components  are  established  and  these 
return  the  segment  to  equilibrium,  (7)  changes  in  the  components  of 
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tension  at  the  end  point  of  the  segments  lead  to  changes  in  the  spring 
forces  in  the  anchors  and  hence  to  changes  in  anchor  locations,  and 
(8)  once  new  equilibrium  positions  of  segments  and  anchors  are  estab¬ 
lished,  the  segments  stop  moving. 

A  new  velocity  magnitude  can  then  be  considered.  If  the  velocity 
increases  or  decreases,  the  static  frictional  resistance  is  recomputed 
and  compared  to  the  holding  capacity  as  described  above. 

Governing  Equations 

Basic  Assumptions .  Several  basic  assumptions  apply  throughout  the 
derivation  of  this  simulation  model.  They  are: 

1.  The  hose  is  uniform  in  shape  and  material  along  the  axis  of 
each  segment  length. 

2.  The  hose  is  infinitely  long  in  comparison  to  its  diameter, 
and  hence  has  no  flexural  stiffness. 

3.  The  hose  is  pressurized  and  maintains  a  circular  cross  section 
without  significant  buckling. 

4.  The  hoseline  slides  over  the  seabed  without  rolling  if  current 
loads  exceed  friction  capacity  of  sea  bottom. 

Governing  Equations.  The  hoseline  system  is  arbitrarily  located  on 
the  seafloor  and  subjected  to  gravity  and  nonuniform  distributed  current 
loads.  Gravity  acts  in  the  direction  of  fc,  which  is  given  by: 

k  =  i/i.  e,  (1) 

ii 

where  the  summation  convention  on  repeated  indices  is  invoked  and  </> ^  is 
the  direction  cosine  of  k.  Hose  behavior  will  be  determined  in  terms  of 
the  location  coordinates  and  the  tension  components  at  a  material  point 
P,  which  is  located  at  an  unstretched  arc  length  Sq  from  origin  along 
the  hose. 
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Unit  vectors  fc  and  n,  which  are  tangential  and  normal  to  the  hoseline 
at  P  respectively,  are  related  to  the  Cartesian  coordinates  by: 


,dX  , 

t  =  l-jf )  e  =  0  e 

\  dS/  a  a  a 


(2) 


where  Greek  subscripts  have  range  1,  2,  where  0q  are  the  direction  cosines 
of  fc,  and  dS  is  the  stretched  differential  length. 


n  =  e  xfc  =  e  -  0  efl 
3  ap  a  p 

where  e^  =  permutation  symbol 


(3) 


=  0  if  a  =  0 

=  -1  if  a  >  3 

=  1  if  a  <  p 

Consequently,  the  Cartesian  base  vectors  as  expressed  in  terms  of  tan¬ 
gential  and  normal  component  fc  and  n  become: 

e  =  0  fc  +  e_  0.  n 

a  a  3a  3 

For  a  hose  section,  the  stretched  differential  length  dS  is  related 
to  its  original  length  dSQ  by: 


ds  ,  , 

-7-  =  1  +  e 

ds 

o 


(4) 


where  e  is  the  strain.  An  elastomer  hose  fortified  with  wire  or  fabric 
reinforcements  usually  stretches  nonlinenrly  under  tension.  The  relation 
between  elongation  and  axial  tension  is  best  described  by  empirical  data. 

In  cases  where  the  required  data  are  not  available,  the  hose  elongation 
is  often  approximated  by: 

e  =  Cj  T  C2  (5) 

where  T  is  axial  tension  magnitude,  while  C^  and  C^  are  material  constants. 

As  the  hose  segment  is  assumed  to  have  no  flexural  or  torsional 
stiffness,  the  tension  acts  in  the  direction  of  unit  tangent  with  magni¬ 
tude  T: 


T  =  Tfc  =  T0  e  =  T  e 
a  a  a  a 
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and  components  of  T  are: 


T  =  T  0 
a  a 


Therefore,  using  Equation  1,  we  obtain  a  first-order  differential  equation 


dX  T 
_  o  _  a 

a  =  dS~  “  T~ 


(6) 


for  in  terms  of  t  and  magnitude  T  given  by: 

T  =  (Tp  •  Tp)1/2 

Replacing  dS  by  dSQ  from  Equation  4,  we  obtain: 
dX  T 

dT  =  a  +  (7) 

o 

Dynamic  Equilibrium.  From  dynamic  equilibrium  of  force  vectors 
acting  on  the  free  body  in  Figure  2: 


+ 

where  T 

fw 


Canceling 


4?]  -  T  +  Ft  dS  +  F„  dS  +  F  dS  +  Fn  dS  +  Fr  dS  =  0 

dS )  L  B  woD  f 

=  tension  at  the  material  point  P  along  the  hose  segment 

=  current  induced  lift  force 

=  current  induced  drag  force 

=  submerged  weight  per  unit  length  of  unstretched  hose 

=  sea  bottom  reaction  normal  to  the  seafloor 
=*  sea  bottom  resistance  parallel  to  the  seafloor 
T  and  dividing  by  dS,  we  obtain: 


dT 

dS 


+ 


fl  +  fb 


+  F 


w 


0 


(8) 
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The  various  force  vectors  are  developed  below. 

Current-Induced  Forces 


Velocity  Vector.  The  current  velocity  acting  on  the  hose  segment 
in  Cartesian  coordinates  is: 

V  =  V  e 

a  a 

which  can  be  expressed  in  normal  and  tangential  components  as: 

V  =  V  + 

n  t 

where  Vfc  =  velocity  tangential  to  segment 
V^  =  velocity  normal  to  segment 

Express  V  in  the  Cartesian  coordinate  system  and  in  the  tangent 
direction  as 


V.  =  V.  e  =  V.  £ 
t  ta  a  t 


where 


Vfc  =  (V  •  £)  i  =  (Vp  0p)  € 


Hence, 


and 


vt  =  Vp 


vt  =  (vpV(0aV 


Therefore, 


Vt«  =  VP  9P  9a 


Similarly,  express  Vn  as: 


V 

n 


V  e 
no  a 


V 

n 


ft 


(9) 


8 


Since 


V  =  V  -  V. 
n  t 


V  =  V  e  =  V  e  -  e 
n  na  a  a  a  to  a 


Therefore, 


V  =  V  -  V.  8.  0  =  Vfl(6  e 

na  a  p  p  a  P  aP 


0o  V 


where  6  .  is  Kronecker  delta 
oS 

6oP  =  °»  a  t  3,  and 


*«9  =  *•  °  "  6 

Also  from  Equation  3: 


Vn  =  (V  •  h)  =  (Vjl  e,)  •  (^  8a  y 
Vp  Eap  0a 

Drag,  (Fd).  Assume  that  the  drag  force  Fp  dS  can  be  expressed  as 
components  in  the  n  and  t  directions: 

Fn  •  dS  =  Fn  •  dS  +  F  _  •  dS 
D  Dn  Dt 

By  the  independence  principle  discussed  in  Reference  4,  tangential 
components  can  be  estimated  in  the  same  manner  as  normal  components 
(Ref  1) 


Fn  =  ■=  p  D  •  Cn  V  •  V 
Dn  2  Dn  n1  n 


(Ha) 


FDt  =  2  P  D  *  CDt  lVtl  '  Vt 


(11b) 


where  cy  and  are  normal  and  tangential  drag  coefficients.  For  a 
hose  segment  of  practical  surface,  C^t  is  much  smaller  than  C^.  From 
Equations  9  and  10,  the  magnitudes  of  | V t |  and  | |  are: 

lVtl  =  <Vta  '  <VS  Va  %  °a)1/2 


|V  I  =  (V  V  )1/2  =  [V-  V  (6.-0  9.)(6  -  8  8  )]1/2 

1  n1  na  na  1  P  n  a3  a  P  an  a  n  1 
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Express  the  drag  force  t-irm  of  Equation  8  in  Cartesian  components: 
fd„  =  (I  o  D  CD„Ki,  •  \/'2  v„« 

FDt  ’  (I  ’  D  CDt)<Vt„  '  Vt,>1/2  Vta 
where  V^a  and  are  given  by  Equations  9  and  10. 

Lift,  (F^) -  dS  acts  in  the  direction  of  -e3>  Assume: 

Fl  -  -  i  P  Ct  D(V0  V„)  S3  (12) 

where  p  =  fluid  density 
D  =  hose  diameter 

CT  =  lift  coefficient 

JLi 

According  to  Reference  7,  the  lift  coefficient  will  remain  positive  but 
will  decrease  as  the  hose  separates  from  the  bottom.  The  hose  will 
repeatedly  be  suspended  in  the  strong  current  momentarily  and  then  fall 
back  to  the  sea  bottom.  This  situation  is,  however,  simplified  in  the 
analysis  by  assuming  the  hose  is  suspended  near  the  sea  bottom,  if  the 
current  is  strong  enough  to  lift  the  hose. 


Gravity  and  Buoyancy  (Submerged  Weight,  F  ) 


Hose  weight  is  specified  as  the  wet  weight  per  unit  unstretched 
length,  w.  The  net  gravity  force  exerted  on  the  hose,  Fw  dSQ,  acts  in 
the  direction  k  and  can  be  expressed  in  terms  of  stretched  length  as: 

/  dS  .  <■ 

F  dS  =  (F  dS)  f-rr2 )  =  dS 

wo  w  \dS  )  1  +  s 

In  Cartesian  components,  we  replace  R  by  Equation  1  to  obtain: 


F  dS 
w  o 


w  dS 
1  +  e 


(4>„  e  +  <P o  e.) 
a  a  j  j 


(13) 


Also,  in  normal  and  tangential  components: 


F  dS 
w  o 


=  dS  (F 


w3 


e_  +  F  n 
3  wn 


+  F„t » 
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where  F  „  =  - 
w3 

w 

L  + 

7  *2 

l 

(14) 

/ 

dS  . 

Fwn  ~  1 

F 

l  w 

dsl 

(•  h  = 

w 

1  +  E 

E«0  \ 

(15) 

t 

dS  , 

W 

1  +  E 

Fwt  =  l 

F 

l  w 

O  1 

dS  J 

Ip  0 
a  a 

(16) 

and  0^  is  given  by  Equation  6  in  terms  of  T  /T. 
Reactive  Forces  From  the  Sea  Bottom 


Reaction  in  e_  Direction,  F_.  As  long  as  the  hose  remains  in 
J  15 

contact  with  the  bottom,  a  reactive  force  F„  dS  will  act  in  the  -e0 

D  J 

direction. 

Considering  equilibrium  in  the  e^  direction,  noting  that  T  has  no 
e3  component  and  using  Equations  12  and  14,  we  obtain: 


'  I  P  D  Cl  V»  V«  -  FB  +  TTH?  =  0 


which  can  be  solved  for  F_ 


B 


W  *3  1 

~  4  p  D  CT  V 
1  +  e  2  K  La 


V 

a 


As  V  increases,  F„  decreases  toward  F„  =  0,  which  occurs  when  the  hose 
o  B  B 

segment  is  lifted  off  the  seabed. 


Frictional  Resistance,  F^.  Assuming  that  the  friction  force  F^  dS 
can  be  taken  as  components  in  the  n  and  £  directions,  and  each  is 
independently  proportional  to  the  bottom  reaction  force,  Fg,  we  have: 

Ff  dS  =  Fg  dS  (vn  n  +  yfc  £)  (17) 

where  y  =  coefficient  of  friction  in  n  directions 
n 

y^  =  coefficient  of  friction  in  £  directions 
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Signs  of  v*n  and  y^.  are  selected  to  oppose  the  tendency  to  move.  Because 
of  the  general  shape  of  hose  segment  under  current  loads,  y^  is  assumed 
to  be  a  skew-symmetric  function  over  the  segment  length 


where  L  =  span  length,  S  =  arc  length  from  the  shore  end,  and  y  is  given. 
Expressing  Equation  17  in  Cartesian  components  by  Equations  2  and  3. 


Ff  dS  =  Ffi  ds(yn  %  +  yt  6po)  0p  ep  (18) 

where  0^  given  by  Equation  6  in  terms  of  T^.  When  the  segment  is  not 
me 
F. 


moving,  the  static  friction  force  F^  is  less  than  the  breakout  force, 


fb 


Ffs  “  I  Fn  ^ 


fb 


s  B 


where  yR  =  coefficient  of  static  friction,  with  |vnl>  | |  <  y 


and  the  signs  of  yn>  yt  taken  the  same  as  Ffg.  Ffg  is  determined  from 
the  static  equilibrium  in  the  normal  direction  as  follows: 


(19) 


If  the  segment  is  not  moving,  the  friction  force  equilibrates  the 
tension  gradient  and  other  external  forces  in  Equation  8  as  long  as  the 
friction  force  is  less  than  the  bottom  resistance  capacity  yF_n. 

D 

Thus,  we  consider  the  n  components  in  Equation  8  to  determine 
static  friction  force  for  comparison  to  the  breakout  force  (holding 
capacity),  ygFB-h.  The  incremental  line  tension  dT/dS  can  be  decomposed 
as  follows: 


dT 

dS 


(  1  )  d(T  £) 

U  +  t)  dso 


t 


+  T 


dt 

dS 

o 
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But 


,?  de  ,  ,T  . 

dt  _  _ a  a  d  |  a]  a 

dS  dS  e«  dS  \T  j  ea 

o  o  o  ' 


_  (1  a  a  dT  I  .  (1  a  a  dT  .  D  a, 

\T  dS  "  T2  dS  )  ea  \T  dS  "  _2  dS  J(0a  t  +  EPa  0p  n) 

'  o  T  o  T  o' 

Therefore, 

dT  =  f  1  \  [:(dT_  dT_\ 

dS  \1  +  t)  z\dS  dS  a  a  a  dS  j 


,dT 

A  (  a 

n  \dS 


o  o 


0P  a  0p  dS 
o  o 


)  ePo 


T  .  T  T  dT 

Since  9a  =  f  ,  0O  0a  =  1>  and  dS-  =  dS^ 

o  o 


T  dT 


■  (rh)  ‘(rss2)-" 


T„  dT  /  T  T 

_P  _ a  fx  _H 

3a  T  dS  \  oti  ’  T  T 
o 


,  _§  _J1  (fi 

EPa  T  dS  l5 


Setting  the  n-components  in  Equation  8  to  zero,  and  using  Equations  20, 
14,  15,  17,  and  19: 


EPo  TP 

dT 

B 

,  T 

X  . 

T  , 
J\) 

w  K  tr 

1  °  r  P 

1  +  E  T 

dS 

o 

\  OTI  T 

T  J 

1  +  e  Pa  T 

+  J  p  D  Cn  (V  V  )l/2  V  6-  s-  +  F.  =0 
2  r  Dn  nri  nil  a  Pa  j  fs 


Solving  for  F^g,  assuming  Ta  is  known,  wo  have: 


T6  *♦.  .1 


eR  *r  T~r~  +  £  P  D  C_  (V  V  )  '  V 
Pa  T  1  +  s  2  v  Dn  nri  nr)  a 


,  dT  ,  T  T 

,  _ 1 _ a  [  r  _a  ji 

1  +  e  dS  an  _  T  T 

ft  ' 
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Breakout  occurs  if 


lFJ  =•  *.h 

in  which  the  sign  of  y  and  F.  =  sign  of  F  .  .  To  evaluate  dT  /dS  in 

S  t)  jlS  T|  O 

Equation  21,  the  actual  computation  uses  central  differences  with  dis¬ 
crete  values  of  at  regularly  spaced  points. 


Dynamic  Equilibrium  in  e^  Direction 


If  static  breakout  occurs,  we  replace  the  static  friction  coeffi¬ 
cient  by  dynamic  friction  coefficient  given  by  Equation  17  and  consider 
equilibrium  in  the  e^,  and  e„  directions.  Set  the  eo  components  in 
Equation  8  to  zero  using  Equations  13,  11,  and  18  along  with: 


dT  _  dT 
dS  dS 


Therefore,  solving  for  dTa/dSo: 

dT  r  Tft  1 

=  -(l  +  e)  ,  ~~ .  +  F„^  +  Fn  +  Fn  ~  (y  ea  +  u^  6„  )  (22) 

dS  v  '  1  +  e  *ct  Dta  Dna  B  T  VKn  Ba  Ht  BcT  v  ' 

o  L  J 

FD„«  -  (I  ”  D  CD„Kn  V„a 

Vna  =  (5oB  '  T  f")VB 
FDt«  ■  (I  »  D  cDt)(vtn  vtn)1/2  Vta 


/dS  .  d(T  e  ) 

[  oi  _  a  a  | 

(  1 

1  ) 

\dS  j  dS  ' 

o 

U  +  ej 

U  +  t) 

T  tr 

V  =  J  J  V 

vta  T  T  B 

Thus,  we  must  solve  Equation  22  as  a  first  order  differential  equation 
along  with  Equation  7,  here  rewritten  as: 

dX  T 

dS  K  }  T 

o 
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for  T  and  X  with 
a  a 

C2 

e  =  C1T 

T  =  (T  T  )1/2 
a  a 

If  the  segment  moves,  we  solve  for  the  new  equilibrium  position  and 
tension  in  which  the  segment  ceases  to  move. 

Numerical  Solutions 


Rewrite  the  governing  Equations  7  and  22  as: 
dX 

_ a  _ 

dS  “  8lct 
o 

dT 

_ a  _  , 

dS  82o  g3a 

o 

Ta 

where  gla  =  (1  +  e)  — 


82o 

83a 


,w  *  > 

■  -(1  +  E>(rT7  +  FD„„  +  FDta) 

’  *16(FB)(l‘n  *B«  +  "t 


The  solution  to  these  nonlinear  equations  can  be  obtained  by 
iteration  on  quasi-linearized  equations.  Taking  Taylor  expansions  of 
g^a>  i  =  1,3,  we  can  use  the  Newton-Raphson  iteration  method  as 
described  in  Reference  8: 


dX 
_ a 

dS 

o 


+  al0BtTa 


dT 

dS~  =  82o  +  83o  +  (a2ag  +  a3a3)(Tp  "  V 
o 
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where  g1#  =  g 


io 


X  ,T 
Tl*  ri 


aiap 


3g . 

&ia 


3T„ 


w  <v 

X^,T^  =  approximated  values 

The  iterative  process  proceeds  as  follows:  assume  values  X  and  T 

at  each  point,  calculate  g ^  and  and  solve  the  linearized  equation 

for  improved  answer  X^  and  T  .  Repeat  this  iterative  process  until 

(X  -  X  )  and  (T  -  T  )  =  0. 
o  o  o  o 


PARAMETRIC  ANALYSIS 

Parameters  Considered  in  the  Test 

A  sensitive  analysis  was  conducted  to  investigate  the  influence  on 
the  hoseline  response  of  the  governing  parameters  considered  in  the 
simulation  model.  The  parameters  tested  include  four  categories:  (1) 
the  axial  rigidity  of  the  hose,  (2)  the  hoseline  geometry,  (3)  the  cur¬ 
rent  velocity,  and  (4)  the  seabed  resistance.  Tests  were  conducted  by 
varying  one  parameter  at  a  time  while  keeping  the  others  constant.  A 
wide  range  of  hose  rigidity  was  tested  for  the  influence  of  rigidity. 
However,  only  a  rigid  hose  was  used  in  the  rest  of  the  numerical  experi 
ments  in  order  to  separate  the  effect  of  other  parameters  from  the  in¬ 
fluence  of  hose  elongation. 

Axial  Rigidity  of  the  Hose 

High-strength  hoses  are  often  made  of  elastomers  fortified  with 
steel  or  synthetic  fiber  reinforcement.  Their  axial  rigidities  depend 
on  the  material  and  the  construction  of  their  reinforcements.  These 
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hoses  usually,  stretch  nonlinearly  under  tension  as  shown  in  Figure  3. 

The  data  presented  in  Figure  3  were  obtained  from  an  elongation  test  of 
a  highly  stretchable  hose  (Ref  9).  The  age  and  the  previous  loading 
history  of  the  hose  may  also  change  the  axial  rigidity  significantly. 
Therefore,  it  is  rather  difficult  to  accurately  estimate  the  rigidity  of 
a  hose.  A  numerical  test  was  conducted  to  examine  the  influence  of  the 
rigidity  on  the  behavior  of  a  hoseline.  In  order  to  study  the  effect  in 
general,  the  rigidity  is  represented  by  a  load  deflection  relation 
approximated  by  Equation  5.  The  influence  of  the  material  coefficients, 
and  C^,  were  tested  separately.  The  test  was  first  conducted  using  a 
hose  of  linear  material.  In  this  case,  is  unity  and  is  directly 

related  to  the  modulus  of  elasticity.  A  wide  range  of  C^,  which  repre¬ 
sents  conduits  varying  from  a  steel  pipe  to  a  highly  stretchable  rubber 
hose,  was  tested.  Figures  4a  and  4b  present  the  lateral  deflection  and 
the  hose  tension,  respectively.  Both  are  nondimens ionalized  by  their 
corresponding  values  of  an  unstretchable  hose,  which  has  a  modulus  of 
elasticity  equal  to  Estee^-  Ehose  anc*  Esteel  rePresent  tlle  n>°dulae  of 
elasticity  of  the  hose  and  steel,  respectively.  S/L  denotes  the  segment 
length-to-span  ratio.  The  results  indicated  that  the  response  of  a  tight 
segment  with  small  S/L  ratio  to  a  strong  current  is  heavily  influenced 
by  the  rigidity  of  the  hose.  This  fact  is  especially  true  when  the 
equivalent  modulus  of  elasticity  is  less  than  one  thousandth  of  Esteej* 

It  is  important  to  note  that  the  rigidity  of  a  typical  submarine  hose  is 
within  this  range.  A  highly  stretchable  hose  in  a  strong  current  elongates 
extensively  and  develops  a  deeper  curved  shape  than  would  a  rigid  hose 
with  the  same  initial  S/L  ratio.  Therefore,  the  tension  load  and  the 
lateral  deflection  of  a  tight  stretchable  segment  in  a  strong  current 
may  increase  as  much  as  300  percent,  if  improper  hose  rigidity  were  used. 

On  the  other  hand,  the  response  of  a  tight  or  loose  segment  in  moderate 
currents  appears  to  be  much  less  sensitive  to  the  variation  of  the  rigidity. 
The  axial  tension  remains  nearly  constant  over  a  wide  range  of  the  rigidity. 
The  reason  can  be  attributed  to  the  effect  of  the  equilibrium  shape  of 
the  segment  on  the  hose  response,  which  will  be  described  in  the  test  of 
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hoseline  geometry.  Figure  5  summarizes  the  effect  of  the  material  coef¬ 
ficient  C^.  The  value  of  reflects  the  stiffening  process  of  the  rein¬ 
forcement  of  a  hose.  For  a  realistic  hose,  is  generally  less  than 
unity.  Again,  the  influence  of  C„  becomes  important  only  when  the  hose 
is  heavily  loaded  beyond  the  extent  that  the  hose  begins  stiffening. 

This  occurs  when  a  tight  hose  is  installed  in  a  strong  current. 

Current  Force 

The  current -induced  force  on  a  hose  section  is  a  function  of  the 
current  velocity,  the  current  direction,  the  hose  diameter,  and  the 
hydrodynamic  coefficients.  This  simulation  model  employs  the  Morison 
equation  to  estimate  the  current  force.  Therefore,  the  current  forces 
can  be  expected  to  be  proportional  to  the  hose  diameter  and  the  hydro- 
dynamic  coefficients.  If  the  curved  shape  of  the  hose  segment  remains 
the  same,  the  force  should  be  proportional  to  the  normal  velocity  squared. 
This  relation  was  verified  by  the  results  presented  in  Figure  6.  The 
hose  tension,  due  to  a  broadside  current,  follows  a  parabolic  function 
closely,  as  long  as  the  current  is  strong  enough  to  overcome  the  seabed 
resistance.  The  influence  of  the  current  direction  on  the  hose  response 
is  more  complicated,  because  the  segment  changes  its  equilibrium  shape 
in  compliance  with  the  current  direction,  as  shown  in  Figure  7.  A  hose 
section  oblique  to  the  current  direction  experiences  a  much  larger  re¬ 
duction  in  form  drag  than  the  increase  in  skin  drag.  As  a  result,  the 
total  current  force  exerted  on  the  section  decreases  significantly. 

Figure  8  summarizes  the  anchor  loads  on  three  segments  of  different  S/L 
ratios  in  various  current  directions.  The  forces  are  nondimens ionalized 
by  the  forces  due  to  broadside  currents.  Generally,  the  anchor  loads 
decrease  in  a  form  of  Cosn8  as  the  current  shifts  away  from  the  direction 
normal  to  the  hose  length.  The  angle  is  defined  in  Figure  9.  The  value 
of  n  decreases  as  the  S/L  ratio  of  the  segment  increases.  Since  the 
hose  is  fairly  stiff  and  stretches  very  little,  the  result  reflects  the 
influence  of  the  current  direction  alone.  Considering  a  segment  of  S/L 
=  1.0,  the  segment  remains  straight  irrespective  of  the  current  direction. 
The  anchor  load  can  be  expected  to  be  proportional  to  the  normal  velocity 
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2  2 

squared  (i.e.,  V.  cos  0),  which  is  confirmed  by  the  results  of  the  case 
of  S/L  =1.01.  On  the  other  hand,  a  loose  segment  aligns  a  large  portion 
of  its  length  with  the  current,  and  therefore  experiences  less  cross 
current.  As  a  result,  the  effect  of  the  current  diminishes. 

Hoseline  Geometry 

The  equilibrium  shape  of  a  hose  segment  determines  the  orientation 
of  the  hose  section  with  respect  to  the  current,  and  therefore  signifi¬ 
cantly  affects  the  total  current  applied  on  a  hoseline.  Intuitively,  a 
loose  segment  tends  to  align  a  larger  portion  of  its  length  with  the 
current  direction,  and  therefore  experiences  a  smaller  current  force 
than  a  tight  segmc  t.  The  curved  shape  of  a  flexible  tension  member 
subject  to  a  uniform  load  can  be  uniquely  characterized  by  the  ratio  of 
the  segment  length,  S,  to  the  span  distance,  L.  Consequently,  the  S/L 
ratios  were  selected  to  represent  the  geometry  of  a  hoseline  segment. 
Figure  10  summarizes  the  influence  of  the  hoseline  geometry  on  the 
hoseline  behavior.  Values  of  S/L  were  generated  by  using  different 
combinations  of  various  span  distances  and  various  segment  lengths.  The 
hoseline  responses  in  1.5-  and  4.0-knot  broadside  currents  were  tested. 
The  results  clearly  demonstrate  that  the  response  of  a  hose  segment  of 
little  stretchability  is  determined  solely  by  the  S/L  ratio  and  is 
linearly  proportional  to  the  segment  length  and  the  span  distance.  In 
general,  the  hoseline  response  remains  roughly  constant  for  a  S/L  ratio 
greater  than  1.2,  and  increases  rapidly  when  the  ratio  decreases  below 
1.2.  Since  an  unstretchable  hose  is  used  in  this  test,  the  lateral 
deflection  is  not  sensitive  to  the  current  velocity  and  the  line  tension 
and  the  anchor  load  are  closely  proportional  to  the  velocity  squared,  as 
anticipated. 

Sea  Bottom  Resistance 

Seabed  resistance  depends  on  the  combination  of  the  unit  weight  of 
conduit,  the  lifting  force,  and  the  soil  properties.  The  soil  reaction 
is  more  complicated  and  less  understood  than  other  parameters  involved. 
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The  mechanism  of  the  soil  reaction  is  approximated  with  a  Coulomb  fric¬ 
tion  force.  All  the  factors  involved  are  included  in  a  simple  friction 
coefficient.  Friction  coefficients  varying  from  0.3  to  0.9,  which  repre¬ 
sent  a  typical  hard  sandy  seabed,  were  tested  to  show  the  influence  on 
the  hose  response.  The  tests  were  repeated  for  conduits  of  unit  weights 
from  3  to  7  pounds  per  linear  foot  in  a  1.5-knot  current.  Unfortunately, 
the  results  shown  in  Figure  11  are  not  conclusive.  Further  research 
efforts  are  required  to  identify  the  influence  of  the  sea  bottom  resis¬ 
tance.  However,  the  typical  submarine  hoses  are  usually  light  and  are 
likely  to  be  lifted  off  the  sea  bottom  by  a  moderate  current.  There¬ 
fore,  the  influence  of  the  bottom  resistance  may  be  negligible. 


STABILITY  OF  NAVY  HOSELINE 

A  three-segment  hoseline  secured  with  four  anchors  equally  spaced 
along  a  straight  line  on  a  flat  horizontal  seabed  was  used  to  show  the 
stability,  in  general,  of  a  Navy  hoseline.  Stability  was  evaluated  in 
terms  of  the  maximum  lateral  deflection,  the  maximum  tension,  and  the 
anchor  loads.  Calculations  were  conducted  for  two  hoses  of  various 
layouts  in  the  current  environment  of  practical  application.  The 
results  were  compiled  in  a  series  of  design  charts. 

Characteristics  of  Navy  Hoses 

Figure  3  demonstrates  the  behavior  of  a  typical  Navy  submarine  fuel 
hose  under  tension.  This  particular  hose  is  made  of  synthetic  rubber 
fortified  with  2-ply  contra  helical  steel  wire  reinforcements.  Figure 
3a  relates  the  axial  elongation  and  the  variation  in  outside  diameter  of 
the  hose  to  the  axial  tension.  The  hose  stretches  linearly  at  low  tension, 
and  becomes  highly  nonlinear  as  the  tension  exceeds  20  kips.  This  fact 
can  be  attributed  to  the  two  distinct  deformation  mechanisms  of  the  hose 
reinforcements  at  different  load  levels.  At  low  tension,  the  helical 
wire  reinforcement  stretches  like  a  regular  helical  spring  without  actually 
deforming  the  wire  material.  Therefore,  the  hose  is  fairly  flexible  and 
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a  linear  relation  between  load  and  deflection  can  be  expected.  As  the 
axial  tension  exceeds  the  magnitude  which  rotates  the  lay  angle  of  the 
reinforcement  wire  to  the  limit,  the  wire  begins  stretching  itself  to 
resist  the  external  loads. ,  The  hose  becomes  much  more  stiff  and  highly 
nonlinear  thereafter.  Another  remarkable  feature  of  the  Navy  hoses, 
along  with  their  high  stretchability,  is  their  extensive  diameter  re¬ 
duction.  Figure  3c  shows  that  the  outside  diameter  of  the  hose  reduces 
almost  the  same  amount  as  the  axial  elongation  in  percentage  of  their 
original  dimension.  Furthermore,  the  results  from  a  simple  calculation 
based  on  the  data  presented  in  Figure  3  indicate  that  the  wall  thickness 
of  the.  hose  changes  very  little  when  the  hose  stretches  (Appendix  B). 
Reduction  in  hose  size  is  due  to  contraction  of  the  inner  diameter.  The 
reduction  in  hose  diameter  has  to  be  properly  accounted  for  in  the  cal¬ 
culation  of  current  loads  experienced  by  the  hose  segments.  The  axial 
stiffness  of  the  hose  is  also  dependent  on  the  pressure  inside  the  hose. 
High  internal  pressure  resists  the  radial  contraction  of  the  hose  rein¬ 
forcement  layer.  This  makes  the  reinforcement  layer  harder  to  stretch, 
as  for  a  helical  spring.  As  a  result,  the  hose  becomes  stiffer  when  the 
pressure  inside  the  hose  increases,  as  shown  in  Figure  3(b).  Therefore, 
the  empirical  load  deflection  curve  of  the  hose  measured  under  working 
internal  pressure  shall  be  used  for  the  final  design. 

Hose  Layouts.  The  calculation  was  repeated  for  two  hoses  of  dif¬ 
ferent  axial  rigidity.  One  is  the  highly  stretchable  Navy  hose  described 
in  the  previous  section.  The  other  is  a  rigid  hose  of  a  very  little 
stretchability.  Hoses  that  were  7.5  and  9.5  inches  in  outside  diameter 
were  used.  The  segment  lengths  of  200,  500,  1,000,  and  2,000  feet  were 
included  in  the  calculation.  The  segment  length  in  combination  with 
different  span  distance  between  anchors  results  in  various  S/L  ratios 
from  1.05,  1.10,  1.15,  and  1.20. 

Environmental  Conditions.  The  hoseline  stability  was  calculated 
for  the  currents  approaching  from  0,  15,  and  30  degrees  off  the  per¬ 
pendicular  to  the  general  direction  of  the  hoseline.  Six  different 
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current  speeds  in  the  0.5-  to  4-knot  range  were  used  for  each  current 
direction.  The  hydrodynamic  coefficients  for  a  submarine  flowline  vary 
significantly  from  case  to  case.  Their  values  depend  on  the  experiment 
setup,  model  scale,  current  condition,  data  acquisition  procedure,  and 
data  reduction  method.  Generally,  these  force  coefficients  are  valid 
only  for  conditions  similar  to  those  for  which  they  were  measured.  The 
inline  drag  coefficients  used  in  the  contemporary  engineering  practice 
vary  from  0.75  to  1.40.  The  transverse  drag  coefficients  vary  from  0.6 
to  1.0.  The  low  values  were  empirical  data  measured  in  37  feet  of  water 
offshore  Honolulu  (Ref  10),  whereas  the  high  values  were  recommended  by 
Det  Norske  Veritas  for  designing  a  submarine  pipeline  system  (Ref  11). 
Both  sets  are  for. steel  or  concrete  coated  steel  pipes.  Data  for 
synthetic  rubber  hoses  are  not  available.  However,  practical  hoses  are 
much  smaller  than  the  thickness  of  the  boundary  layer  of  the  water  cur¬ 
rent,  and  are  therefore  fully  submerged  in  the  boundary  layer.  Con¬ 
sequently,  the  actual  current  velocity  experienced  by  the  hoses  is 
smaller  than  that  of  the  free  stream.  Using  the  lower  force  coeffi¬ 
cients  in  combination  with  the  free  stream  velocity  seems  to  be  more 
realistic.  Navy  hoseline  is  a  relocatable  system  intended  for  worldwide 
application.  The  system  may  be  installed  on  various  types  of  sea  bot¬ 
toms.  Field  data  show  that  the  typical  Coulomb  friction  coefficient  for 
a  steel  or  concrete  pipe  on  a  hard  sandy  seabed  varies  from  0.3  to  0.9 
(Ref  5).  A  low  friction  coefficient  of  0.3,  which  will  result  in  higher 
loads,  was  used  in  the  calculation.  The  hose  particulars  and  the  en¬ 
vironmental  conditions  are  summarized  .in  Table  1. 

Results 

A  total  of  eight  design  charts  (Tables  2  through  5)  were  prepared 
for  various  hose  rigidities,  hose  particulars,  and  current  conditions. 
The  results  for  the  rigid  hose  (Tables  2  and  3)  were  further  reduced 
into  a  parametric  model.  Equations  23  through  26,  using  the  findings 
obtained  from  the  parametric  analysis  as  a  guideline.  The  parametric 
model  is  compared  with  the  complete  simulation  model  in  Figure  12.  The 
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solid  lines  indicate  two  models  with  identical  results.  The  parametric 
model  provides  a  good  approximation  for  the  conditions  within  the  limits 
indicated: 
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deflection  (ft) 

maximum  hose  tension  (lb) 

current  load  at  end  anchors  (lb) 

current  load  at  intermediate  anchors  (lb) 

span  distance  between  anchors  (ft) 

segment  length  between  anchors  (ft) 

normal  component  of  current  velocity  (knots) 

outside  diameter  of  the  hose  (in.) 


Limits: 


5  <  D  <  10  inches  200  <  S  <  2,000  feet 

15V  <4  knots  1.05  <  S/L  <  1.20 

n  ' 

Engineering  Applications 

ACP  Hose.  The  stability  of  the  Advanced  Collapsible  Pipe  hose  can 
be  evaluated  directly  by  using  the  design  charts.  It  is  assumed  that  a 
6-inch  ACP  hose  is  to  be  installed  in  a  uniform  2-knot  cross  current. 
The  outside  diameter  of  the  hose  is  approximately  7.5  inches.  It  is 
further  assumed  that  the  system  is  to  be  installed  at  an  initial  (or 
unstretched)  segment  length-to-span  ratio  of  1.05.  The  anchor  capacity 
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required  to  hold  a  hose  segment  of  various  lengths  against  a  2-knot 
broadside  current  can  be  obtained  from  Table  3  under  the  column  labelled 
with  '0  degree'  and  'Anchor  load.'  Columns  (1)  and  (2)  give  the  current 
loads  at  the  end  and  the  intermediate  anchors,  respectively.  For  example, 
the  minimum  holding  capacity  of  the  anchor  at  the  end  of  a  1,000-foot 
segment  can  be  obtained  from  the  row  labelled  with  S/L  =  1.05,  Length  = 
1000,  and  V  =  2.0.  The  results  are  3.2  and  4.2  kips  for  the  end  and  the 
intermediate  anchor,  respectively.  As  a  matter  of  fact,  there  is  more 
than  one  way  to  anchor  a  hoseline.  Final  selection  requires  a  trade-off 
analysis  taking  into  consideration  of  the  strength  of  the  hose,  the 
lateral  deflection,  the  holding  capacity  of  the  anchors,  the  installa¬ 
tion  requirements,  and  the  operation  costs,  etc. 

Other  Hoses.  The  stability  of  other  hoses,  whose  axial  rigidities 
are  substantially  different  than  those  of  the  ACP  hose,  may  be  deter¬ 
mined  by  iteration  using  the  design  charts  for  the  rigid  hose  along  with 
the  empirical  load  deflection  relation  of  the  hose  under  consideration. 

The  iteration  procedure  includes  four  simple  steps:  (1)  determine  the 
hose  tension  from  the  proper  design  chart  for  the  rigid  hose  using  the 
initial  (unstretched)  S/L  ratio  and  other  given  parameters,  (2)  deter¬ 
mine  the  elongation,  e,  associated  with  that  tension  from  the  empirical 
load  deflection  relation,  and  (3)  calculate  the  new  S/L  ratio  by  multi¬ 
plying  the  old  S/L  by  (1+e).  Then,  (4)  use  the  new  S/L  ratio  and  repeat 
steps  1  to  3  until  the  tension  load  comes  within  desired  tolerance.  For 
example,  the  stability  of  the  ACP  hose  demonstrated  in  the  previous  para¬ 
graph  can  also  be  determined  using  the  design  charts  for  rigid  hoses. 

Step  1:  Entering  Table  2  with  S/L  =  1.05,  V  =  2.0,  S  =  1,000,  and  0=0 
degree,  the  tension  load  is  4.5  kips.  Step  2:  Interpolating  from  Figure 
3,  the  elongation  of  the  hose  under  4.5  kips  tension  is  7  percent.  Step 
3:  The  new  S/L  ratio  is  therefore  equal  to  1.12.  Repeating  steps  1  to  3 
with  the  new  S/L  ratio,  the  second  iteration  gives  a  tension  of  2.8  kips 
and  a  elongation  of  5  percent.  The  third  iteration  gives  a  tension  of 
3.1  kips  and  an  elongation  of  5  percent,  which  are  within  the  accuracy 
of  the  design  charts.  Therefore,  the  equilibrium  S/L  ratio  of  the  seg¬ 
ment  is  1.05  X  (1+0.05)  =  1.10.  The  loads  at  the  end  and  the  intermediate 
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anchors  are  determined  from  Table  1,  with  S/L  =  1.10  and  S  =  1,050  feet, 
as  3.2  and  4.2  kips,  respectively.  The  results  are  identical  with  those 
obtained  directly  from  the  design  charts  for  ACP  hoses.  Furthermore, 
the  stability  of  a  rigid  hose  can  also  be  determined  from  the  parametric 
model,  Equations  23  through  26.  Equation  24  is  used  in  place  of  design 
charts  for  estimating  the  tension  loads,  and  Equations  23,  25,  and  26 
are  used  to  determined  the  lateral  deflection  and  the  anchor  loads.  The 
iteration  procedure  remains  the  same. 


CONCLUSIONS 

1.  The  lateral  stability  of  a  flexible  submarine  hoseline  on  the  sea¬ 
floor  in  a  slowly  varying  current  environment  may  be  properly  simulated 
based  on  the  Morison  equation  and  a  nonlinear  cable  theory.  Stability 
can  be  evaluated  in  terms  of  anchor  load,  axial  tension,  and  maximum 
lateral  deflection  of  the  hoseline. 

2.  The  on-bottom  behavior  of  the  hoseline  is  most  influenced  by  the 
following  factors:  the  equilibrium  curved  shape  of  the  hose  segment, 
the  size  of  hose,  the  current  velocity,  and  the  axial  rigidity  of  the 
hose.  The  influence  of  the  sea  bottom  resistance  is  negligible  from  a 
design  point  of  view. 

3.  The  equilibrium  curved  shape  of  a  hose  segment,  which  may  be  repre¬ 
sented  by  the  segment  length-to-span  ratio,  S/L,  is  the  most  dominant 
factor  on  the  response  of  a  hoseline.  The  influence  of  S/L  on  the 
responses  is  illustrated  in  Figure  10.  The  results  show  that  a  tight 
hose  segment  with  a  small  S/L  ratio  w.ilJ  experience  large  current  loads 
with  small  lateral  deflection,  and  that  a  loose  hose  with  a  large  S/L 
ratio  will  experience  a  smaller  current  load,  but  with  a  much  larger 
lateral  deflection.  A  proper  S/L  ratio  has  to  be  selected  according  to 
the  capacity  of  anchors,  the  strength  of  the  hose,  and  the  operational 
requirements.  A  ratio  between  .1.05  and  1 . 10  is  recommended  for 
practical  applications  based  on  the  results  of  a  parametric  analysis. 
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4.  At  a  constant  equilibrium  S/L  ratio,  the  load  is  proportional  to  the 
hose  size,  the  hydrodynamic  coefficients,  and  the  span  distance,  and  is 
closely  proportional  to  the  normal  current  velocity  squared. 

5.  The  axial  rigidity  influenced  the  elongation  of  a  hose  segment,  and 
hence  the  equilibrium  S/L  ratio.  The  influence  is  significant  for  a 
tight  stretchable  hose  segment  deployed  in  a  fast  current,  and  is  negli¬ 
gible  for  other  cases  such  as  a  loose  hose  of  little  to  moderate  stretch- 
ability  deployed  in  a  slow  current. 

6.  The  optimum  segment  geometry  of  a  simple  hoseline  may  be  determined 
by  iteration  using  the  design  charts  provided  in  Tables  2  to  5  or  the 
parametric  model  described  by  Equations  23  to  26.  The  iteration 
procedures  are  demonstrated  in  the  section  on  Engineering  Applications. 
The  simulation  program,  NCELHOSE,  is  recommended  for  the  analysis  of  a 
more  complicated  hoseline.  The  source  code  and  assistance  to  use  the 
program  are  available  at  NCEL. 
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Table  1.  Range  of  Parameters  Used  to  Generate  Design  Charts 


Parameter 

Range 

Hose  Diameter 

Outer 

7.5  and  9.5  inches 

Inner 

6.0  and  8.0  inches 

Hoseline  Geometry 

Length 

200,  500,  1,000  and  2,000  feet 

S/L 

1.05,  1.10,  1.15,  and  1.20 

Current  Velocity 

Magnitude 

0.5,  1.0,  1.5,  2.0,  3.0,  and  4.0  knots 

Direction 

0,  15,  and  30  degrees 

Hydrodynamic  Force 
Coefficients 

& 

0.75 

0.60 

Ct 

0.1 

Friction  Coefficients 

Static 

0.3 

Dynamic 

0.25 

t 
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Table  2.  Response  of  rigid  hoseline  to  steady  currents  (Do  =  7.5  inches,  RIGID 

Di  =  6.0  inches.  Wet  =  5  lb/ft,  =  0.75,  CT  =  0.6,  =  0.3)  - 
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Table  2.  Continued 
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Table  3.  Response  of  rigid  hoseline  to  steady  currents  (Do  -  9.5  inches, 
Di  =  8.0  inches,  Wet  =  5  lb/ft,  Cn  =  0.75,  C.  =  0.6,  =  0.3) 
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Profile 


Gravity  Vector  g  k  SHORE 


X2  Direction 


Vector 


Figure  1.  Definition  sketch. 


FdtdS  FRdS 


Figure  2.  Vectors  on  a  typical  hose  segment. 


(a) 

o 


(b)  (c) 


Figure  3.  Behavior  of  a  highly  stretchable  hose  under  axial  tension.* 
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Nondimensionalized  Tension  Nondimensionalized  Deflection 


Figure 


Figure 


5.  Influence  of  material  constant  C2  on  the  response  of  a  hose 
in  a  4 -knot  current. 


Current  velocity,  V  (Knots) 


6.  Influence  of  current  velocity  on  the  hose  tension,  S/L  =  1.2. 
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Figure  7.  Geometry  of  the  hose  segment  in  currents  of  various  directions. 


Current  direction  (Degrees) 


Figure  8.  Influence  of  current  direction 
on' the  anchor  load,  V  =  4 
knots. 


Hose 


Figure  9.  Definition  of  current 
direction. 
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Anchor  loads/Length  (lb/ft)  Max.  axial  teiision/Length  (lb/ft)  Max.  deflection/len 


Segment  length-to-span  ratio  Segment  length-to-span  ratio  , 


Figure  10.  Influence  of  segment  length-to-span  ratio,  S/L,  on  the 
response  of  a  hoseline. 
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Numerical  model  Numerical  model 

10  15  20  25  0.0  0.1  0.2  0.3  0.4 


Parametric  model  (Eq.  23) 


Parametric  model  (Eq.  25) 


Tin 

L 

□  / 

□o/te 

Parametric  model  (Eq.  24) 


Parametric  model  (Eq.  26) 


Figure  12.  Parametric  model  versus  numerical  model. 


Appendix  A 


NUMERICAL  SOLUTION  TECHNIQUES 


A.l  NEWTON-RAPHSON  QUASI -LINEARIZATION 

Given  a  coupled  set  of  nonlinear  first-order  differential 
equations,  it  is  possible  to  develop  a  convergence  acceleration 
procedure  (Ref  A-l)  of  successive  iteration  upon  quasi-linear  equations 
adopting  methods  proposed  by  Bellman  (Ref  A-2)  and  by  Klaba  (Ref  A-3). 
Assume  a  set  of  2N  nonlinear  first-order  differential  equations: 


=  f^s.y^  i,j  =  1,2,...,2N 


(A.1.1) 


with  N  boundary  condition  at  the  starting  end  s  =  0 


gm(yj)  =0  m  =  1,2, ... ,N 


(A. 1.2a) 


and  N  boundary  conditions  at  the  terminating  end  S  =  1, 


h  (yT)  =0  m  =  1,2,. . . ,N 
m  wj  ’  5  ’ 


(A. 1.2b) 


where  S  is  the  independent  variable  (e.g.,  hose  arc  length),  y^  are  the 
2N  dependent  variables  (e.g.,  tension  components  and  location  coordinates), 
f^  (s,yj)  are  nonlinear  functions  of  y^  (e.g.,  hydrodynamic  forcing  func¬ 
tions),  and  g  (y?)  and  h  (y^)  are  nonlinear  combinations  of  y?  at  s  =  0, 
and  of  y^  at  s  =  1,  respectively. 

Let  y^  denote  a  trial  solution  vector  in  the  neighborhood  of  the 
true  solution  vector  y ^ .  The  y?1  and  y^ '  are  corresponding  boundary 
values  of  yl  at  s  =  0  and  s  =  1,  respectively.  The  nonlinear  function 
f^(s,yj),  gm(y?),  and  fm(yl)  can  be  written  as  truncated  Taylor  series 
expansions  about  y.',  y?'  and  y^'  as: 


fi(s,yj)  =  fi<s>yp  +  Jik(yk  "  yk}  i,j,k  =  1,2,...,2N  (A.1.3) 


2  (y°)  =  e  (y°')  +  J°  (y°  -  y°') 

®nryj'  ®nryj  1  mk^yk  yk  ' 


m  =  1,2, ... ,N 


h  (yl)  =  h  (y^1 )  +  J*.  (yf  -  y?' )  m  -  1,2, . . . ,N 

m wj  mVJj  mk  Jk  Jk  >  >  > 

where  the  summation  conversion  from  1  to  2N  on  repeated  indices  has  been 
adopted,  is  the  square  Jacobian  matrix  of  order  2Nx2N  for  the 
gradient  of  the  nonlinear  forcing  function: 


3fi(s,yi) 


yj  yj 


(A. 1.4) 


and  and  J^  are  the  rectangular  Jacobian  matrices  of  order  Nx2N  for 
the  gradients  of  the  nonlinear  boundary  conditions: 


J°  =  J° 
mk 


3g  (y?) 
°m  J  i 


(A. 1.5) 


jl  =  jlk 
mk 


J  J  J 


(A. 1.6) 


Upon  substitution  of  Equations  A.1.3  and  A. 1.4  into  Equation  A. 1.1 
and  A. 1.2,  the  boundary-value  problem  can  thus  be  written  as: 


ar  =  aik(s)  yk +  bi(s) 


(A. 1.7) 


with  boundary  conditions 


0  0  ,  ,0  _  _ 
c  .  y  .  +  d  =0 
mk  mk  m 


(A. 1.8a) 


A-2 


1 

Cmk 


where 


1  .  ,1  _ 
y  .  +  d  = 
mk  m 

0 

(A. 1.8b) 

lik 

Jik 

(A. 1.9a) 

b.  = 
i 

fi(s,yi) 

•  Jik  *k 

(A. 1.9b) 

,0 

mk 

3  O 
7T 

(A. 1.9c) 

d°  = 
m 

,0  0, 

Jmk  yk 

(A. 1. 9d) 

'mk 

A 

mk 

(A. 1.9e) 

d1  = 
m 

h  (yl)  - 
m  Jj 

T1  1. 

Jmk  ^k 

(A. 1.9f) 

Equation  A. 1.7  with  boundary  conditions  A. 1.8  and  coefficients  defined 
by  A. 1.9  constitutes  a  linear  boundary  value  problem  for  an  improved 
solution  y^  in  terms  of  functions  for  the  previous  trial  h^.  Further 
improved  solutions  are  obtained  by  successive  iteration  in  Equations 
A.  1.5  to  A.  1.9  with  yj.,  y0',  and  y1,  replaced  by  the  y,,  y^,  y*  gen¬ 
erated  by  the  previous  iteration.  The  iteration  process  continues  until 
the  difference  between  y^  and  yl  is  less  than  a  stipulated  accuracy. 


A. 2  DECOMPOSITION  OF  LINEAR  BOUNDARY-VALUE  PROBLEM 

A  linear  two-point  boundary-value  problem  such  as  that  posed  by 
Equations  A. 1.7  and  A. 1.8  can  be  solved  by  first  decomposing  the  problem 
into  a  set  of  initial-value  problems  and  then  recombining  solutions  to 
each  initial-value  problem  (Ref  A-4,  A-5,  A-6).  The  advantage  of  using 
this  method  for  solving  a  two-point  boundary-value  problem  rather  than 
a  different  method  is  that  large  sets  of  matrix  equation  coefficients 
need  not  be  generated,  stored  in  the  computer  mem.'ry,  and  solved  simul¬ 
taneously.  Only  a  small  number  of  coefficients  at  the  starting  and 
terminating  points  need  to  be  considered.  According  to  Ince  (Ref  A-7), 
the  solution  to  each  one  of  a  linear  set  of  2N  first-order  differential 
equations  can  be  considered  as  a  linear  combination  of  the  solutions  of 
N+l  initial-value  problems,  hereafter  called  partial  solutions. 


A-3 


Assume  the  solutions  to  Equation  A. 1.7  can  be  written  as: 


0  ,  n 
yi  =  4  •. 


(A. 2.1) 


where  e^,  n  =  1,2,...,  N  are  undetermined  parameters  and  z^  and  za  are 
partial  solutions.  The  partial  solution  z?  is  the  particular  solution 
to: 


i  0  , 

ds  aik  Zk  +  bi 


(A. 2. 2a) 


subject  to  actual  "initial"  conditions  at  s  =  0 


cO  ,  zO,  +  dO  =0 
mk  k  m 


(A. 2. 2b) 


and  fictitious  "initial"  conditions  at  s  =  0 


cO  .  z0.°  +  d0°  =  0 

mk  k  m 


(A. 2. 2c) 


The  partial  solutions  z^,  n  =  1,2,...,N  are  homogeneous  solutions  to 


dz . 

i  _  n 

ds  aik  Zk 


n  =  1,2, . . . ,N 


(A. 2. 3a) 


subject  to  actual  "initial"  conditions  at  s  =  0 


cO  ,  zO  =  0 

mk  k 


n  =  1,2, . . .,N 


(A. 2. 3b) 


and  fictitious  "initial"  conditions  at  s  =  0 


cO  .  zO?  +  d0n  =  0  n  =  1,2, ... ,N 

mk  k  m  >  >  > 


(A. 2. 3c) 


The  choice  of  fictitious  initial  conditions  coefficients  cO  ,  ,  dO  , 

n  n  mk  m 

d0m  is  such  that  d0m  are  linearly  independent  vectors  and  cO^  is  a 

rectangular  Nx2N  matrix  of  coefficients  which  allows  an  inverse  of  the 

assembled  2Nx2N  square  matrix  C ^  of  initial  value  coefficients: 


j,k  =  1,2,. ..,2N 
m  =  1,2,. ..,N 


(A. 2. 4) 


A-4 


to  obtain  solutions  for  zO 


0 

k 


and  zO?1 
k 


as: 


m 


(A.2.5) 


m 


(A. 2. 6) 


The  dOn  are  usually  taken  as  Kronecker  delta  functions  and  dOn  as  a  null 
vector. 

Having  defined  N+l  linearly  independent  initial-value  problems, 
each  of  which  satisfies  the  actual  boundary  conditions  at  s  =  0,  one 
integrates  each  problem  to  the  terminating  point,  s  =  1.  The  partial 
solutions  obtained  at  the  terminating  point  are  then  used  to  determine 
the  appropriate  parameters,  en>  in  Equation  A. 2.1  for  the  linear  com¬ 
bination  of  partial  solutions.  The  boundary  conditions  expressed  by 
Equation  A. 1.8b  at  s  =  1  can  be  written  in  terms  of  partial  solutions 
as: 


cl»k(zlk  +  +  dln,  =  0  (A'2'7) 

0  n 

where  zl^  and  zl^  are  the  terminal  values  of  the  partial  solutions.  The 
product  cl^zl^  is  a  square  NxN  matrix  and  thus, 

en  =  -Kk  *«S'lK  +  cl»*  21°)  <a-2-8> 

With  en  determined,  a  final  initial-value  integration  of  Equation  A. 1.7 
can  be  performed  with  initial  values: 


y0k 


zO?  +  e  zof 
k  n  k 


(A. 2. 9) 


A-5 
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Appendix  B 

REDUCTION  IN  DIAMETER  OF  AN  OPDS  HOSE  UNDER  TENSION 


The  Navy  Offshore  Petrolem  Discharge  System  (OPDS)  relies  on  a 
6-inch  (inner  diameter)  rubber  hose  for  fuel  transfer.  The  hose  is  made 
of  multilayer  elastomeric  material  fortified  with  two-ply  contra-helical 
wire  reinforcements  to  withstand  the  external  tension  loads.  The  hose 
stretches  to  almost  45  percent  of  its  original  length  at  fracture 
(Figure  3  in  main  text).  In  the  meantime,  the  outside  diameter  of  the 
hose  reduces  40  percent.  This  significant  size  reduction  tends  to  shut 
down  the  fuel  flow  and  totally  disable  the  hose line.  Research  for  hose 
construction  techniques  to  control  the  hose  neck-down  is  currently 
underway,  which  will  be  addressed  in  the  final  project  documentation  of 
the  Advanced  Collapsible  Pipe  Program.  This  analysis  is  to  examine  the 
change  in  the  inner  diameter  of  the  hose  under  tension.  Elastomer  is 
generally  treated  as  an  incompressible  material.  Its  total  volume  does 
not  change  when  subjected  to  external  forces  (Ref  B-l).  Referring  to 
the  definition  sketch.  Figure  B-l,  the  total  volume  of  the  hose  wall,  V, 
can  be  expressed  as: 


v  =  n  D  t  Jt 


in  which  D  =  (D  +  D,)/2 
o  i  ' 


where 


D 

t 

H 


outer  diameter 
inner  diameter 
mean  diameter 

wall  thickness  of  the  hose 
total  length  of  the  hose 


Since  the  elastomer  is  assumed  incompressible, 

AV  AD  ,  At  ,  Afc  _  n 
-V  ='d  +  “  +  T  -  0 


(B-l) 


where  A  is  the  differential. 

Replacing  the  mean  diameter  D  with  Dq  -  t  and  assuming  that  Dq  »  t: 

D  =  D  -  t 
o 

=  Do  •  -i > 
o 

AD  _  o  _  n  ,  _tv  f _ o_  At 

D  D  (1  -  t/D  )  U  D  ;  \  D  "  D 

o'  '  o  o  o  o 

B-l 


Substitution  .AD/D  into  Equation  B-l  and  neglecting  the  second  order  term 
At/Do> 

At  _  [A*  ,  M  ,  t  ,  ADol  r 

-  -  -  K  +  (1  +  -)  — I  (B-2) 

L  o  o  J 

Therefore,  the  variation  of  the  wall  thickness  under  tension  can  be 
expressed  in  terms  of  the  elongation  and  the  outer  diameter  reduction  of 
the  hose.  Table  B-l  summarizes  the  results  of  Equation  B-2  using  the 
empirical  data  shown  in  Figure  3  in  the  main  text.  The  far  right  column 
shows  that  the  variation  in  the  wall  thickness  is  less  than  6  percent, 
even  though  the  outer  diameter  of  the  hose  reduces  more  than  40 
percent.  This  result  implies  that  a  substantial  reduction  occurs  in  the 
inner  diameter. 
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Table  B-l.  The  Diameter  reduction  of  the  Navy  Hose  Under  Tension. 


Tension 

(kips) 

O.D. 

(in) 

dl/t 

dD  /D 
o  o 

l+t/D0 

dt/t 

1 

0.0 

7.8 

- 

~ 

- 

10.0 

6.9 

0.0734 

-0.0548 

1.1176 

-0.0122 

20.0 

6.1 

0.0940 

-0.1159 

1.1229 

0.0262 

30.0 

5.6 

0.0469 

-0.0820 

1.1441 

0.0469 

40.0 

5.3 

0.0261 

-0.0536 

1.1643 

0.0363 

50.0 

5.0 

0.0109 

-0.0566 

1.1799 

0.0559 

60.0 

4.9 

0.0144 

-0.0200 

1.2013 

0.0096 

70.0 

4.8 

0.0106 

-0.0204 

1.2074 

0.0140 

80.0 

4.7 

0.0070 

-0.0208 

1.2147 

0.0183 

90.0 

4.7 

0.0174 

0.0000 

1.2233 

-0.0174 
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